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Abstract 

The exact solution of the ID interacting mixed Bose-Fermi gas is used to calculate ground-state 
properties both for finite systems and in the thermodynamic limit. The quasimomentum distri- 
bution, ground-state energy and generalized velocities are obtained as functions of the interaction 
strength both for polarized and non-polarized fermions. We do not observe any demixing instability 
of the system for repulsive interactions. 
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The cooling and trapping of quantum Fermi gases of ultracold atoms poses a number of 
additional challenges to those faced for bosons [1]. A key point is that no more than one 
identical fermion can occupy a single state due to the Pauli exclusion principle. However, a 
mixed gas of fermions and bosons provides an effective means of cooling single-component 
fermions by thermal collisions with evaporatively cooled bosons [2-5], providing an avenue for 
investigating many-body quantum effects in degenerate Fermi gases. Another development 
is the use of Feshbach resonance, in which the energy of a bound state of two colliding 
atoms is magnetic field tuned to vary the scattering strength from — oo to oo, allowing the 
investigation of the crossover from BCS superfluidity to Bose-Einstein condensation [6-8]. 
These achievements in realizing quantum Fermi gases of ultracold atoms may also provide 
insights into other areas of physics, such as ultracold superstrings [9]. 

Particular theoretical and experimental interest has been paid to the Fermi gas confined 
in ID geometry [10-14]. Recent attention has turned to the ID model of mixed bosons 
and polarized fermions [15-18], revealing quantum effects such as interaction-driven phase 
separation [17-19], bright solitons in degenerate Bose- Fermi mixtures [20] and Luttinger 
liquid behaviour [15, 16]. In this model, only s-wave scattering for boson-boson and boson- 
fermion interactions is considered, with the boson-fermion collisions minimizing the Pauli- 
blocking effect in momentum space. 

The theoretical studies of the mixed Bose- Fermi gas have focused on the case of bosons 
and single-component fermions [15-18]. Our aim here is to investigate the ID model of a 
Bose-Fermi mixture on a line of length L with periodic boundary conditions, with Hamil- 
tonian H = — + Yli<j 9i,j$( x i — x j)i where among the N particles there are 

Nf = + Af fermions and iVj, bosons. The pairwise (^-interaction has strength g it j. The 
crucial observation from a theoretical point of view is that if all particles have equal masses 
and if the interaction strength between all particles is the same, the above model is exactly 
solvable by Bethe Ansatz (BA) [21]. Although this restricts the applicability of the model, 
experiments with isotopes of atoms are expected to meet the BA conditions [18]. Then 
g iy j = 2h 2 /(ma), where c = 2/a is the inverse ID scattering length of the confined particles. 
For convenience of notation, the energy is measured in units of h 2 /(2m), such that 



The model contains two special limiting cases: (i) the ID interacting Fermi gas with ar- 
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bitrary polarization [22-25], and (ii) the ID interacting Bose gas [26]. The mixed model was 
recently discussed for the case of single component (fully polarized) fermions [18]. Here we 
consider the more general case of two-component fermions with arbitrary polarization. We 
present analytical and numerical results for the ground state energy, quasimomentum den- 
sity profile, and the spin and charge velocities. In the weak coupling limit, these quantities 
reveal the typical signatures of pure quantum gases with an additional weak interaction due 
to the mixture. In the other extreme, that is for strong repulsive interactions, the ground 
state properties resemble those of a single-component non-interacting Fermi gas. 

The BA equations (BE), determining the quantum numbers of the A-particle system, 
are given by [21] 

rpi - Aj + id 
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where d = c/2 and j, k = 1, . . . , Aj, + A]-, £ = 1, . . . , N, m,n = 1, . . . , A&. Without loss 
of generality, we take Aj < N^. Here the set of N + Aj + 2A^ many quantum numbers 
is divided into three subsets: {pe}i,...,N, {Ak}i,...,N r +N b , {Ai}i,...,iv 6 - It turns out that the 
energy eigenvalues E are given by the members of the first set alone, namely E = J2f =1 pj- 

We begin with a finite system, for which the Bethe roots are found analytically in the 
weak coupling limit, thereby yielding the ground state energy. The second step is to carry 
out the thermodynamic limit (TL). In this limit, the Bethe roots are distributed smoothly 
over a certain interval of the real axis, giving rise to continuous densities. Integral equations 
for these densities have been obtained previously [21]. On the one hand, these equations 
allow a comparison between the analytical results for weak coupling and finite systems with 
the weak coupling expansion in the TL. On the other hand, generalized velocities can be 
calculated within this framework. In the weak coupling limit, these reduce to the spin and 
charge velocities of a pure Fermi gas and the charge velocity of a pure bosonic system. 

In considering the ground state for weak interaction, it is convenient to distinguish be- 
tween unpaired pf ] (j = 1, . . . , N t - A T ), paired pf ] (j = 1, . . . , A T - 1, A T + 2, . . . , 2A T ) 
and bosonic p^ (j = 1, . . . , N b + 2) quasimomenta [27]. Expanding Eqs. (2)-(4) to 0(c), 
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one obtains pf = tt(-A^ - 1 + 2j)/L + 5^, j = 1, . . . , - iV T )/2; pj u) = 7r(2iV T - N t - 1 + 

(p) 



■l-Ni + 2j + ) + 8f>±y/c/L. 



2j)/L + 5j u) , j = {N i - iV T )/2 + 1, . . . , N i - 7V T and pf = vr( 
Here j + = j if j odd and j + = j — 1 if j even. The deviations 5j from ^ u ' p - ) are linear in c, 
with 
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Here Pj P (j u ^ denotes the quasimomenta in the free particle limit as given above. The sums 



in (5), (6) over paired momenta p^ count each pair only once. The bosonic momenta p 
collapse at the origin for c = 0. Their c-dependence is given by 
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According to (7), the p^ are the roots of Hermite polynomials of degree N b + 2 [28, 29]. 

From the above equations we obtain the ground state energy, E = E<®+E<U, to leading 
order in c. The energy of the free particles, E°, is given by the corresponding expression for a 
free Fermi gas with Nf fermions [25]. The first correction in c is = + E^p + 2j^N b Nf, 
where E^ (E^) is the linear order in c for a pure Bose (Fermi) gas of N b bosons [29] (Nf 
fermions [25]). The last term is due to the interactions between bosons and fermions. Thus 
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E = 



2tt 
~L 



-JV T (JV? - 1)+ 



^(N i - JV T )(-1 + Nf + N 2 f - N f N,)j 
+J + l N b(N b - 1) + N b N f ^j . 
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This expression is valid for both weak repulsive and attractive interaction. Especially, the 
TL is well defined in the weakly attractive case, as opposed to the strongly attractive case 
[23]. Let us now carry out the TL, i.e., N a , L — > oo where the densities n a = N a /L are 
held constant, with a f, b. In this limit, Eqs. (5) and (6) give the distribution of 
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quasimomenta per unit length, 

^ P) =l + ^{^ + ?)' A<lPl<C W 
1 c / A C 2B\ 

Pp ^~n ~ 2^ U 2 "P 2 + C 2 -p 2 ~ ~W) ' 

I^cBJk <\p\<C (10) 

Pb(p)~ (4c5/7r-p 2 ) 1/2 , \p\ < 2y/tffa, (11) 

where A = vriV T /L, 5 = nN b /L, C = tt(N 1 - N^/L. 

The quasimomentum distribution function of the bosonic momenta (11) for c > is given 
by the semi-circular law behaviour as for a pure bosonic system [26, 28, 29]. For c < 0, 
the p^ in (7) are imaginary, so that (11) with imaginary k yields the dark-soliton like 
distributions of the bosonic quasimomenta on the imaginary axis. These encode the binding 
energy of the bosons, a quantity potentially accessible by experiments. In the distribution 
of the fermionic quasimomenta, divergences occur near the cutoffs. The quasimomentum 
distribution functions calculated from the BE (2)- (4) are compared with the approximations 
(9)-(ll) in Fig. 1. 

In order to study the effect of arbitrary interaction in the TL, we use the linear integral 
equations derived by Lai and Yang [21] which determine the densities of the quasimomenta. 
In the following, we restrict ourselves to two limiting cases, namely TVf = N± and TVj = 0. 
In the latter case, the fermions do not interact among themselves due to the Pauli principle, 
so that the interaction potential in (1) is only effective between fermions and bosons. The 
ground-state energy for this case has been calculated recently [18]. In the other limiting 
case, iVf = Ni, all the fermions interact with each other and with the bosons. Fig. 2 shows 
the ground-state energy per unit length E/L for different densities rib and rif obtained 
from numerical solution of the integral equations. Also shown is a comparison between the 
analytic result (8) and the TL. 

Up to this point, we have focussed on the density p = p u + p p + p b of the quasimomenta 
p( u 'P' b )_ j n an analogous fashion, one may introduce densities a, r of the roots Aj, A m . Using 
the dressed energy formalism [23], one can calculate the corresponding dressed energies 
€,<j),ijj, which are the energies necessary to add or remove a root to or from the seas of 
Pi,Aj,A m . The dressed energies give rise to generalized velocities, which determine the 
asymptotics of correlation functions [30]. In the case of pure bosons, there is only one 
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species of BE numbers, associated with one dressed energy function, yielding the charge (or 
sound) velocity vjp [23, 26]. For pure fermions, the two sets of BE numbers are linked with 
the charge and spin velocities vi{] . As has been proven by Haldane [31], these velocities 
coincide with those calculated in the harmonic- fluid (or bosonization) approach [32]. 

The situation is less clear for the mixture of bosons and fermions. In the BA approach, 
this corresponds to two seas of BE numbers pe, Aj, and correspondingly two dressed energies 
and two velocities, which we call vi b ' f \ For c = we have = 27rn/ and vf^ = 0. The 
dependence of vf'^ on c for different values of nf,n b is shown in Fig. 3. The lowest c-order 
of vi b ^ is obtained from the Bethe-Ansatz as = 2n b [c/n b — (l/27r)(c/n b ) 3 ^ 2 ] which is 
the velocity of a non-interacting Bose gas for small c [26]. 

We now compare our results to those of the hydrodynamic (HD) approach [15]. In this 
framework, the pure gases correspond to harmonic oscillators, such that the coupling between 
them leads to new normal modes v±, with v±(c = 0) = v^' b \c = 0). In the weak interaction 

_ 2 

limit, V- = 2y/au, v+ = 2nnf + For the sake of comparison, the latter result is 

also shown in Fig. 4. In the strongly repulsive limit, a demixing instability is predicted 
from the HD [15] and mean-field approaches [17]. We do not observe any instability for 
repulsive interaction for the integrable model, in agreement with [18, 21]. The reason for 
the discrepancy is that to our present understanding, the HD approach for the mixture, 
especially the calculation of normal modes, is a low-energy theory. That is, it is expected 
to yield reliable results for small interaction strengths |c| <C 1. Investigation of v¥' b ^ within 
the BA approach in the strongly repulsive limit yields 



where n — rif + n b . 

It remains a very interesting question for future research if the normal modes of the 
field theory coincide with the generalized velocities from the BA for mixed Bose-Fermi 
systems. Furthermore, it should be carefully investigated whether or not the HD approach 
is applicable in the strongly interacting regime. 

The generalized velocities can also be computed from the BA for interacting fermions 
(n-j- = rijj. In this case, the three dressed energies e,(p,ip give rise to three velocities, which 
we call vi{] and vf\ As shown in Fig. 4, we observe that in the weak interaction limit, 
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VcJ = 7m/(l ± c/(n/7T 2 )), which corresponds to the charge and spin velocities of a pure 
fermionic system [12]. In the strong interaction limit, vi b ^ = vP = 0, whereas vP = 2-rrn. 
This again indicates the fermionic nature of the system in the strongly interacting limit. 

In conclusion, we have investigated ground state properties of the ID interacting Bose- 
Fermi model from its exact BA solution. We obtained results for the distribution of quasi- 
momenta and the ground state energy both for weak attractive and repulsive interactions. 
We computed the generalized velocities and compared them to the normal modes obtained 
from the HD approach. In contrast with other approaches [15, 17, 33], we do not observe 
any instability or demixing in the system. 
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FIG. 1: The quasimomentum distribution functions for a finite system (crosses) (iVf = 15, iVi = 
31, Nb = 9) and in the TL (full lines in the insets) for (top) weakly repulsive and (bottom) weakly at- 
tractive interactions. The insets a),c) show the distributions near the noninteracting cases (dashed 
lines) 1/27T for unpaired fermions and 1/tt for paired fermions. The inset b) depicts the bosonic 
quasimomenta for c > 0; for c < 0, their real part is zero (thick line). 
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FIG. 2: Ground state energy per unit length E/L in the TL as a function of the interaction 
strength for densities = n±,n = n^ + nj = 1, with = 1,0.9,0.7, . . . ,0.1,0 from bottom to top. 
The dotted line is the asymptotic result E/L — > 7r 2 /3n 3 for c — > oo, which is equal to the energy 
of free fermions with density n. The inset compares the analytic (8) (straight dashed lines) with 
numerical TL-results (solid lines) for weak repulsive interaction. 
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FIG. 3: Velocities vi (upper set of curves) and Vc (lower set of curves) for polarized fermions 
with rif = 0.9, 0.7, 0.5, 0.3 (from top to bottom) interacting with bosons n& = 0.1, 0.3, 0.5, 0.7 (from 
bottom to top). The insets show the weakly repulsive limit, together with analytical results from 
BA and HD. 
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FIG. 4: Velocities v y c n /{im } ) (upper curve) and v s /(irrif) (lower curve) for a mixed Bose-Fermi 
gas with rif = rib = 0.25 plotted against c/rif. The upper inset shows the whole c-range. The 
second inset shows for weak couplings. 
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